Goldstone models in D + 1 dimensions, D = 3, 4, 5, supporting 
stable and zero topological charge solutions 
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Abstract 

We study finite energy static solutions to a global symmetry breaking Goldstone 
\ model described by an isovector scalar field in D + 1 spacetime dimensions. Both 

CN ■ topologically stable multisolitons with arbitrary winding numbers, and zero topological 

charge soliton-antisoliton solutions are constructed numerically in D = 3, 4, 5. We 
have explored the types of symmetries the systems should be subjected to, for there 
to exist multisoliton and soliton-antisoliton pairs in D = 3, 4, 5, 6. These findings are 
■ underpinned by constructing numerical solutions in the D < 5 examples. Subject to 

axial symmetry, only multisolitons of all topological charges exist in even D, and in 
odd D, only zero and unit topological charge solutions exist. Subjecting the system 
to weaker than axial symmetries, results in the existence of all the possibilities in all 
dimensions. Our findings apply also to finite 'energy' solutions to Yang-Mills and 
r~| \ Yang-Mills-Higgs systems, and in principle also sigma models. 

<N 

> ■ 1 Introduction 

Very early in the history of field theory solitons, interest in the existence of zero topological 
charge solutions arose. In the case of the Yang-Mills (YM) instantons pQ in D = 4 Euclidean 
space, which are selfdual, this raised the question of existence of non-selfdual [21 H] solu- 
tions, while even earlier this question was investigated [5] in the case of magnetic monopoles 
of the YM-Higgs (YMH) model in D = 3 [6] . More recently, concrete numerical construc- 
tions of monopole-antimonopole solutions [3, El [9] to the YMH model in D = 3, instanton- 
antiinstanton solutions 10 to the YM model in D = 4, as well as soliton-antisoliton solu- 
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5h I tions [TT] to a Goldstone model in D = 3, were given. 

The potential relevance of field theory soliton-antisolitons in higher dimensions rests in 
the fact that they describe non BPS field configurations that may be useful in the description 
of brane-antibrane configurations. Non BPS configurations are relevant for example in the 
context of string junctions in N = 4 super- Yang-Mills [12] ■ Such solutions can be the zero 
topological charge counterparts of higher dimensional instantons [13] and of monopoles p3] , 
or of the solitons of the symmetry breaking Goldstone type models [15] arising as the gauge 
decoupling limits of higher dimensional monopole models [H]. These Goldstone models 
have not found any physical applications to date, but as prototype systems modeling higher 
dimensional monopoles without the burden of gauge degrees of freedom, they can be useful 
for example in providing backgrounds on which Dirac equations [16] in all dimensions can 
be solved, or possibly also for gravitating monopoles. Here, they will prove very useful in 
studying zero topological charge solutions in higher dimensions. 
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Zero topological charge solutions to such a Goldstone model in D = 3 were recently given 
in [11]. The model in [11] is the gauge decoupling limit of the YMH model descending from 
the p = 2 member of the YM hierarchy introduced in [13] and is the simplest example. In the 
present paper we will extend this study to the Goldstone model descending from the p = 3 
member of the YM hierarchy. In contrast to the p = 2 Goldstone model which supports 
finite energy solitons only in D = 3, the p = 3 Goldstone model enables us to study solutions 
in dimensions D = 3,4,5, allowed by the Derrick scaling rule. This is very important for 
our purposes here as will be explained below. As such, the p = 3 Goldstone model will serve 
as a vehicle for us to investigate zero topological charge solutions in the simplest possible 
technical setting. 

The main objective of this work is to find out: subject to what symmetries and for 
what boundary conditions do such solutions exist? We have presented several numerically 
constructed solutions in dimensions D < 5, by way of underpinning our findings. While 
our symmetry considerations cover the dimensions 3 < D < 6, the concrete numerical 
constructions are limited to D = 3, 4, 5 in the p = 3 Goldstone model, covering both even 
and odd dimensions, allowing us to make a classification of the said conditions. Our study 
addresses the question as to what are the requisite ingredients in the construction of zero 
topological charge solitons in higher dimensions, highlighting the distinction between even 
and odd dimensions in this respect. We find that such solutions can be accommodated by 
imposing the requisite boundary conditions for systems subject to the appropriate symmetry, 
in all dimensions. Stated most succinctly, subject to axial symmetry only multisolitons of 
arbitrary charges exist in even D, while in odd D zero and unit topological charge solutions 
can exist. By imposing less stringent symmetries than axial, all possible types of solutions 
can be constructed in any dimension. 

The symmetries considered are at one extreme axial, namely spherical symmetry in a 
R D_1 subspace of R D , and at the other azimuthal, namely rotational symmetry in a R 2 
subspace of R D . In between, we have explored the imposition of all intermediate cases, 
namely imposition of rotational symmetry in all the other subspaces R ra of R D . In addition, 
we have considered imposition of multi- azimuthal symmetries on all R 2 subspaces of R D . 

Concerning the numerical constructions, our reason for limiting to the p = 3 Goldstone 
model, and to D < 5, is that otherwise it would be necessary to carry out numerical inte- 
grations in more than 2 dimensions, which is beyond the scope of this work. 

In Section 2 we have introduced the models to be employed, along with the topological 
charge densities providing the lower bounds on the energies. Section 3 is concerned with the 
imposition of symmetries, i.e. stating the axial, azimuthal, intermediate and multi- azimuthal 
Ansatze. In particular, the energy density functionals of the model, for the dimensions in 
which numerical solutions will be constructed, are subjected to the spherical, axial and 
bi-azimuthal symmetries. Subjecting the corresponding topological charge densities to sym- 
metries is carried out in Section 4. Section 5 contains all the numerical results, which verify 
the assertions presented in the previous Section, 4, concerning the symmetry properties that 
zero topological charge solutions must have. In Section 6 we summarise our results and 
extend the discussion of symmetries to beyond the particular simple models employed here. 
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2 The model and the topological charge 



The symmetry breaking models in D = 3, 4 and 5 spatial dimensions, to which we refer as 
Goldstone models, are described by a scalar isovector field a , a = 1,2,3, a = 1,2,3,4 and 
a = 1, 2, 3, 4, 5 in each dimension, respectively. 

There is such a hierarchy of models [15] that arise from the gauge decoupled limit of 
the D— dimensional SO(D) gauged Higgs (YMH) model descended from the p-th member of 
the Yang-Mills (YM) hierarchy on Rd x S Ap ~ D . Here we have chosen the simplest of these 
that can accommodate D = 3,4,5, while satisfying the Derrick scaling requirement for the 
existence of finite energy solutions. In the present case this is the YMH model that descends 
from the p = 3-rd member of the YM hierarchy. Our Goldstone model here is the gauge 
decoupled limit of this YMH system. Using the notation 

with the brackets [/xja..] implying total antisymmetrisation, the static energy density is 

£ {P =3) = A (v 2 - |0T) 6 + Ax (v 2 - |0T) 4 M 2 + A 2 W - l0 b D 2 \Cf + ^ \C C P \ 2 ■ (!) 

All the dimensionless constants Aq, Ai, A2 and A3 must be positive if the relevant topological 
lower bounds in each dimension are to be valid. The model (CQ) is ad hoc rather than a 
dimensionally descended, only inasfar as the numerical values of these dimensionless coupling 
constants, which are otherwise fixed by the descent mechanism, are constrained only to be 
positive. Of course, any of these constants can also vanish, provided that the absence of the 
corresponding term in ([1]) does not violate the Derrick scaling requirement. 

The most important feature of models such as (PQ) is that the order parameter field 4> a 
is a relic of a Higgs field and has the same dimensions {L~ l ) as a connection, and the finite 
energy conditions allow the symmetry breaking boundary conditions 

lim |0 a | = , lim |0 a | = 77 , (2) 

R being the radial coordinate in R D , this resulting in monopole-Wke asymptotics for our 
solitons. 

The presence of the symmetry breaking potential in ([!]), multiplying Ao, has no quanti- 
tative effect on the solutions, so it will be ignored henceforth. 

In the next Section, where symmetries will be imposed, we will concentrate only on the 
terms 

|tf , |tf| 2 » Itfpf (3) 

and will delay the incorporation of the factors (rf — |0 6 | 2 ) 2 and (rf — |0 fe | 2 ) 4 till the Section 
on numerics, since imposition of symmetries on these last terms is achieved rather trivially. 

The topological charge density bounding the energy density functional from below can be 
stated simply in terms of Bogomol'nyi inequalities, separately for each dimension D = 3, 4 
and 5. 
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In D = 3 the inequality 
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Each of the three topological charge densities g 3 , £> 4 and g 5 is a total divergence, which 
we denote as g 3 = d^Sl^ , f? 4 = d^Vt^ and g 5 = d^Q^, respectively, the surface integrals 
of yielding the topological charge in each dimension D = 3,4,5. In this paper we will 
refer to the densities fijf^ as topological currents. Now these topological charges are simply 
numerical multiples of the respective winding number densities 



(o) 
Qd 



aia 2 ...a D iai ia 2 xa D _ o (D) 



which are the surface integrals of the winding number currents 



to 



in 
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an 



(9) 



(10) 



1 Note that the D = 3 model employed here is slightly different from that in [IT], the latter being the 
gauge decoupled version of the p = 2 YMH model, in contrast to the gauge decoupled version of the p = 3 
YMH model here. 
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The topological charges ^3,^4 and q$, which are the volume integrals of the densities gs, 
£4 and £5 defined in (JSJ), (JTj) and (jSj) respectively, are in turn equal to the surface integrals 
of the topological currents 
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e^ T e abcde L 2 - |0| 2 j <f> a €<f> c p <Pt€ ■ 



(11) 



It is now obvious, in light of the asymptotic boundary value in (j2J), that (73,(74 and q$ are 
multiples of the winding numbers, namely the surface integrals of the currents fflUj) . with the 
numbers « and =, respectively 

3 Imposition of symmetries 

This section is divided into four subsections, in each of which the three building blocks in (J3]) 
will be subjected to spherical, axial, azimuthal, and in 4 dimensions only, the bi-azimuthal 
symmetries respectively. We shall also state the tri-azimuthal Ansatz, only in 6 dimensions, 
but will not display the building blocks subject to it because in 6 dimensions the Derrick 
scaling requires an octic term beyond these, whose numerical pursuit is beyond the scope 
of this work. The symmetry breaking selfinteraction potential (rf — |0 a | 2 ) 6 will be ignored, 
instead the boundary condition (|2]) it would enforce will be imposed directly. 

Imposition of symmetry is the first step in the construction of zero topological charge 
solutions, leading to the second step of selecting the requisite boundary conditions to achieve 
this aim. In this section, we impose the symmetries on the energy density functional ([I]) 
whose second order equations will be integrated numerically in Section 5, deferring the task 
of imposing symmetries on the topological charge densities and their currents (ITU1) to 
the next section, 4. There, the most important task of selecting the requisite boundary 
conditions will be made. 

Before stating the Ansatze pertaining to the various symmetries to be imposed on the 
scalar field <p a describing the model (HI) , we introduce the coordinates to be employed in 
each case. Next to spherical symmetry, the strongest symmetry that we will impose is the 
axial symmetry, sometimes described also as cylindric symmetry. This involves imposition 
of rotational symmetry in a R 0-1 subspace of the full space R D . The weakest symmetry 
is the azimuthal one, which involves imposition of rotational symmetry in a R 2 subspace 
of R D . Then there are all the intermediate symmetries involving imposition of rotational 
symmetry in a R n subspace of R°, with D — 2 > n > 3. As we restrict to D — 5, the only 
relevant values of n are n = 3 and 4. In addition, we will employ multi- azimuthal symmetries, 
each one of its constituent azimuthal symmetries being imposed on distinct planes in R°. 
Since we will restrict to D = 6, our attention will be restricted to the bi-azimuthal and 
tri-azimuthal cases only. The coordinates are parametrised as follows. 
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Axial coordinates : 

In this case we label the coordinate on R D as follows, 

Xfj, = (xi, x D ) , i = l,2,..,D-l , \xi\ 2 =r 2 , R 2 = r 2 + x 2 D , (12) 

so that 

r = Rsm9i , xd — Rcos9\ , (13) 

where 9i is the leading polar angle in each dimension, parametrised by the spherical polar 
angles (#1, 9 2 , #0-3, 9 D _ 2 , <p), <P being the azimuthal angle (with < <p < 2n, < 9i < ir) 
Our definition of axial symmetry amounts to spherical symmetry in the D — 1 dimensional 
subspace, as for example in [T7j. 

Azimuthal coordinates : 

Imposing azimuthal symmetry in the Xi = (x\, Xz) plane, leaves the dependence of the fields 
on the coordinates xj = (#3, £4, --,xd) unrestricted. In practice however, we will restrict to 
the D = 4 case only for reasons explained in Section 4 below. The labeling we will employ 
is 

x ll =(x i ,x I ) , i = 1,2 , 1 = 3,4 , \xi\ 2 = p 2 , R 2 = p 2 + \x I \ 2 , (14) 
so that 

p = R sin 61 sin 9 2 , £3 = R sin 9\ cos #2 , £4 = -R cos 9\ , (15) 

or 

p = r sin ^2 , £3 = z = r cos ^2 , X4 = t . (16) 
Intermediate coordinates : 

In D = 5 the only intermediate possibility is n = 3, and we label the coordinate as 

Xfj, = (xi, x 4 , x 5 ) = (xi, s, t) , i = 1, 2, 3 , \xi\ 2 = r 2 , i? 2 = r 2 + s 2 + t 2 , (17) 
so that 

r = i? sin 6*i sin 9 2 , s = R sin 6 1 ! cos 6*2 , t = Rcos9i, (18) 

in an angular parametrisation (# 1; 9 2 ,93,ip), with polar angles ranging from to it, cp being 
the azimuthal angle ranging from to 2n. The notation 

Xi = (sin 9 3 cos ip, sin 9 3 sin cos 6*3) (19) 

will be employed below. 

In D = 6 both n = 3 and n = 4 are possible, but the second leads to a 4 dimensional 
effective system which is superfluous for our purposes here. Hence we restrict to n = 3 and 
label the coordinate as 

x fi = (Xi,x 5 ,x s ) = (Xi,s,t) , i = l,2,3,4 , \ Xi \ 2 = r 2 , R 2 = r 2 + s 2 + t 2 , (20) 
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so that (r, s, t) are parametrised exactly as in ([15]) in an angular parametrisation 
# 2 , #3, #4, 93). The notation 

Xi = (sin #3 sin # 4 cos if, sin #3 sin 64 sin sin #3 cos 64, cos #3) (21) 

being employed for this case below. 

Bi-azimuthal coordinates : 

In this case we will restrict our attention to D = 4 and D = 5. (Bi-azimuthal symmetry 
in D = 6 would lead to four dimensional residual subsystems, which are superfluous for 
our purposes here.) In the first case we will subject the components of the energy density 
functional ([3]) to the symmetry implied by the Ansatz, while in the second, we will only state 
the Ansatz since no solutions will be constructed subsequently in that case. 

In D = 4 we impose a second azimuthal symmetry in ffT41) . in the xj = (x3,X4), plane, 
denoting the radial variable in the (x, y) and (z, t) planes with p = yj x 2 + y 2 = \/\xi\ 2 and 
a = y/z 2 + t 2 = yj\xi\ 2 . In this case, we will parametrise R 4 as 

Xi = (Rsinip) Xi=pxi , Xi = (cos^i, sin<£>i) (22) 
xj = (Rcosip) xj = oxj , xj = (cosy? 2 , sin<^ 2 ) 

where R 2 = \xi\ 2 + |x/| 2 = |x^| 2 , with < ip < § , < <pi < 2ir and < y? 2 < 2tt. While the 
two angles (ipx, <£> 2 ) are azimuthal angles, the angle 9 here is not a polar angle as its range is 
one half of that of a polar angle. We shall refer to such angles as semi-polar henceforth. 
In D — 5, R 5 is parametrised as 

Xi = (Rsm9smip) Xi = r X{ , Xi = (cosy^, sin</?i) (23) 

xi = (R sin 9 cos ip) xj = s xi , xi = (cos ^2, siny2 2 ) 

x 5 = Rcos9 = t (24) 

where R 2 = r 2 + s 2 + 1 2 , and < 9 < tt and < if) < f . In ([23])-(f24]), 9 is a polar angle and 
ip a semi-polar angle. We denote polar angles by 9 and semi-polar angles by ip henceforth. 
All azimuthal angles are likewise denoted by ip. 

Tri-azimuthal coordinates : 

Here we restrict our attention to D = 6 only for reasons explained already. Extending the 
labeling Q22D of R 4 to that of R 6 , with p = ^Jx\ + x\ = \Z\x~\ 2 , i x = 1, 2, a = sjx\ + x\ = 
y/\xi 2 \ 2 , 12 = 3, 4, and r = a/x 2 + x\ = \/\x iz \ 2 , i 3 = 5, 6, by 

x h = (R sin ipx sin ip 2 ) x h =px h , x^ — (cos^i,sin^i) 

x i2 = (R sin ipi cos ip 2 ) = f7 ^i 2 j £j 2 = (cos</> 2 , siny? 2 ) (25) 

x i3 = (i^cos^i) x i3 =TXi 3 , x i3 = (cosy? 3 ,siny2 3 ) 

where R 2 = (xjj 2 + \xi 2 \ 2 + \xi 3 \ 2 = |a^| 2 , with < ipi < | , < ip 2 < f , and with the three 
azimuthal angles < ipi < 2tc, < ip 2 < 27r and < (^3 < 27r. 
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3.1 Spherical symmetry 

The spherically symmetric Ansatz for the scalar field cj) a in D dimensions is 



x" 



^ = r]Q(R)x a , x a = — 

tx 



(26) 



resulting in the reduced building blocks 
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where we have used the notation Q R = |2 . 



3.2 Axial symmetry 

The axially symmetric Ansatz for the scalar field 4> a 
the index a — 1, 2, ...D — 1 is 



in D > 4 dimensions, with 



^ a = r]H(r,x D )x a , (j) D = r]G{r,x D ) , £ a 



(28) 



j4 ,0 3 ), A = 1,2, is tantamount to imposing azimuthal 



using the labeling ([12]) of the coordinates. 

There is a very important exception in the D = 3 case of ( l28l) . where imposition of 
axial symmetry on the field a = 
symmetry. The axially symmetric Ansatz in D = 3 is 

(f) A = r\ H(r,x 3 ) n A , (J) 3 = r]G(r,x 3 ) , n A 



cosny? 
sin rap 



(29) 



n = 1,2, 3, ... being the azimuthal vortex number. 
The result of substituting (1281) in (|3J) is 
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(tf 2 + G 2 + # 2 + G 2 ,) + (D - 2) , (30) 

2(H [r G D] ) 2 + (D - 2) 2 (tf 2 + G 2 + if 2 + G 2 D ) + (D - 3) ^ 
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6(# [r G D] ) 2 + 



+(D-3)^y 3(^ + G? + flS + G2,) + (D-4)^y | 



where we have used the notation H r = Hp = J^, and H[ r GD] = {H t Gd — HqG 7 ). 

The spherically symmetric limit ([27)1 of (I3"0j) follows immediately from (TlBl) . by the re- 
placements 

if (r, x D ) = Q(R) sin 6' jD _ 2 , G(r, x D ) = Q(R) cos 6 D - 2 , 
with i? 2 = r 2 + a;|,, and using 

a a a , cos ®D~2 o o fl o sin6 l z)_2 „ 

d r = sin 6»£)_2 o R H — de D _ 2 , d D = cos 6> D _ 2 »r — — d dD _ 2 . 

3.3 Azimuthal symmetry 

This subsection is concerned with the imposition of azimuthal symmetry in a D dimensional 
system, resulting in a D — 1 dimensional residual subsystem. As such, it does not lead to 
a boundary value problem which can be tackled numerically in a practical way. It should 
thus be viewed as a first step towards the imposition of bi-azimuthal symmetry in the 
D = 2 + 2 = 4 case presented in the next subsection. 

Imposing azimuthal symmetry in the Xi = (21,22) subspace (plane) of x^ = (xi,xi), 
I = 3, 4, D, and labeling the scalar field as <fi a = (<fi A , <f) A '), A = 1, 2 and A' = 3, 4, £>, 
the components A are restricted by the Ansatz 

<p A = h(p,Xj) n A , n A = (cos rap, sin n<£>) , p 2 = |xj| 2 = x 2 + y 2 , (31) 

while the D — 2 components </> A ' = <p A ' (p, xj) retain their dependence on the D — 2 coordinates 

X/. 

The result of enforcing the Ansatz (I3TI) is most compactly expressed by employing the 
coordinate xu = (%i,p), and by labeling the residual field as x a — (x A iX Dl ) = (0 A ;^)> 
with the new index running over a = A', D — 1. In this notation we have 
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In the case of interest here, namely for D = 4, x M = (xj,x/), with i — 1,2 — x,y and 
7 = 3, 4 = 2, i, the azimuthally symmetric Ansatz fl3~Tj) now becomes 



A = %,^,t)n A , A ' = x A '(x 7 ,p) 



f(p,z,t) 



(33) 



resulting in the residual 3 dimensional system with coordinates % = (z, t, p) being given by 
fl32|) with = X 3 = h. 
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The axially symmetric limit (1301) of (I32I) follows immediately from (fl5l) . by the replace- 
ments 

h(z, t, p) = H(r, t) sin 0i , f(z,t,p) = H(r,t)cos9 1 , g(z,t, p) = G(r,t) , 
with r 2 = p 2 + z 2 , and using 

3.4 Intermediate symmetries 

The symmetries to be considered here are rotational symmetry in the R 3 subspace of R 5 (for 
D = 5), and the R 4 subspace of R a (for D = 6). (Rotational symmetry in the R 3 subspace 
of R 6 would be superfluous since that would lead to a four dimensional effective system.) 
For D = 5, the intermediate symmetric Ansatz for the field a = (0 a , 4 , 5 ) is 

r = r ] h(r,s,t)x a , <f> 4 = V g(r,s,t) , ^ = V f(r, s, t) (34) 

using the notation of (IT7j) - (fTUl) . 

The intermediate symmetric Ansatz for the field <p a = (<f) a , 5 , 6 ) for D = 6 is 

<f> a = rih(r,s,t)x a , (f> 5 = rig(r, a ,t) , 6 = V f(r, s, t) (35) 

which looks formally identical to f[3~4"|) . but now the coordinates being read from (|20p - (|21|) . 

In both cases the system reduces to a three dimensional effective subsystems, for which 
numerical constructions are outside the scope of this work. Hence we do not present the 
result of symmetry imposition on the energy density functionals (j3J). 

3.5 Bi-azimuthal symmetry 

Our considerations in this subsection cover two cases, namely to state the bi-azimuthal 
Ansatze in D = 4 and D = 5. The residual subsystem in each case is two dimensional and 
three dimensional, respectively. In the first case we will construct the solutions numerically, 
so the Ansatz will be imposed on the energy density functional, while in the second we will 
limit ourselves to stating the Ansatz. 

Bi-azimuthal symmetry in D = 4 : 

In the D = 4 case, using the notation (|22|) for the coordinates and using the same notation 
(I3ip as in subsection 3.3, a = ((p A ,(j) A ), the bi-azimuthally symmetric Ansatz is 

4> A = r]h(p,a)n A , n A = (cosni</?i, sinni^i) , (36) 
4> A ' = VOiPi 17 )^' > n 2 — ( cos ^2^2, sin n 2 ^2) , 

where ni and n 2 are the respective vorticities in the two planes. 

In fact the Ansatz (|36|) results in the first stage from the imposition of azimuthal symme- 
try (I3"T|) in D = 4, with the residual fields \ a — (0 A > h), and then imposing a second stage 



10 



of azimuthal symmetry on the triplet x a - Concerning the imposition of the second stage of 
azimuthal symmetry, we point out that the densities (1321) resulting from the first stage do not 
exhibit a global SO(D-l) invariance, although the original densities ([3]) are invariant under 
a global .SO (.0)0. We have verified that the second stage results a consistent reduction, even 
though the reduced system after the first stage did not possess a global invariance. 

Imposition of bi-azimuthal symmetry enables a 2 dimensional boundary value problem, 
to be tackled numerically in the next section, so we list the resulting densities ([3]) 
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where we have used the notation h. 



= dh h _ - 
RsinO, a = 



dh 



and g [p g a] = (h p g a - g p h a ) as in 
In terms of the coordinates p = Rsin.8, a = Rcos8 defined by (I22p . the spherically 
symmetric limit fl30|) of fl37j) follows immediately from by the replacements 



h(p,a) = Q(R) sm9 , g(p, a) = Q(R) cos 9 



and using 
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This limit will be exploited in the numerical constructions. 

Bi-azimuthal symmetry in D = 5 : 

Here, the residual system being three dimensional, we simply state the Ansatz 



A5 



r)h(r,s,t)nf , nf = (cosni</?i, sinni^i) , 
r)g(r,s,t)ri2 , = (cos 71,2^2, sin n 2 (p2) 



(38) 
(39) 
(40) 



in the notation of 



3.6 Tri-azimuthal symmetry 

As noted at the start of this section, we shall simply state the Ansatz here, for 6 dimensions 
only, without imposing the symmetry on the energy density building blocks ([3]). Then in the 
next section we will use this to calculate the topological charge of the putative solutions in 

2 This is in contrast to that of a YM system, where the local gauge group does, under azimuthal symmetry 
imposition, reduce to an effective YM-Higgs system exhibiting a broken local gauge invariance [10] . 
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6 dimensions, which are not constructed numerically here. The tri-azimuthally symmetric 
Ansatz for a = (</A, Aa , As ), A 1 = 1,2, A 2 = 3,4, A 3 = 5, 6 

4> Al = h(p, a, t) nf 1 , nf 1 = (cosni^i, sinrii^i) 

4> A2 = g(p,&, T ) n 2 2 j n;f 2 = (cosn 2 <^ 2 , sinn 2 y? 2 ) (41) 

As = f(p,v,T)ri2 3 , n^ 3 = (cosri 3 y2 3 ,sinri 3 y2 3 ) 

where ni, n 2 and n 3 are the respective vorticities in the three planes (xi,x 2 ), (x 3 ,x 4 ) and 
(x 5 ,x 6 ). 

4 Topological charges and boundary values 

In this Section, we present in detail the topological charges resulting from the various types 
of boundary values of the scalar field. This is relevant when subjecting the fields to axial, 
azimuthal intermediate bi-azimuthal and tri-azimuthal symmetries in turn. Under each 
(symmetry) heading, we will calculate the topological charges in all dimensions D for which 
the residual subsystem is at most three dimensional. This will cover the generic cases, 
all further examples being superfluous. Subject to axial symmetry, we consider the cases 
D = 3, 4, 5, 6. Subject to azimuthal symmetry, we cover only D = 4. Subject to intermediate 
symmetry, we take the cases D = 5, 6. For configurations with bi-azimuthal symmetry, we 
cover D = 4,5. Subject to tri-azimuthal symmetry, we cover the only possible case D = 6. 

As explained at the end of Section 2, it is sufficient to calculate the winding numbers 
since the topological charges are simply numerical multiples of the latter. Up to angular 
volume normalisations A^^, these are the surface integrals of the currents (jTOl) . hence what 
we need to calculate are the asymptotic values of the quantities i^wf to enable us to 
evaluate the surface integrals 

/2 D - 1 ^(0 D _ 2 ,0 D _ 3 ,..,^ 1 ,^), (42) 

R=oo 

being the unit vector, and <KI(9d-2, @d-3, V 9 ) t ne angular volume element, in H D . 
Here, we will evaluate the angular integrals (j4"2l) . (a) subject to axial symmetry for 
D = 3, 4, 5, 6, (6) subject to azimuthal symmetry for D = 4, and (c) subject to bi-azimuthal 
symmetry for D = 4. 

4.1 Axial symmetry 

In the case of axially symmetric fields, we will impose the following asymptotic boundary 
values on the functions H(r,xr>) and G(t,xd) defined in (|28p for D > 4, and in (T2T?|) for 
D = 3 

lim H(r,XE>) = shimf^ (43) 

R^oo 

lim G(r, i D ) = cosm^i , m = l,2,3,... 

R— >oo 



id 
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The topological charges of the axially symmetric models in D = 3, 4, 5 are defined by 
the integrals (H2|) . divided by the angular volumes Qd-i = 2tt,2tt 2 ,^- in each of these 
dimensions respectively, by the (volume) integrals 

q D = =D\ f H D - 2 (G R H 9l - H R G 6l ) dRd6 l , 

S'£)-l J 

The surface integrals ( 1421) can be evaluated analytically. In the axially symmetric cases at 
hand, where the corresponding volume integrals are two dimensional, these become contour 
integrals in the positive half plane r[0, oo), xd(— oo, +oo) by virtue of Stokes' theorem. 
Now the line integral along the xd axis does not contribute since analiticity requires that 
H(6i = 0) = H(9i = 7r) = 0, so the only contribution comes from the infinite semicircle, 
thus reducing fj42l) to the following one dimensional angular integrals 

J axial = D\n D J H D ~ 2 (GH 01 -HG 9l ) \ R=oo d9, , (44) 

with the exception of the D = 3 case where axial symmetry coincides with azimuthal sym- 
metry, when 

h = 2\2ixn J H (GH e - HG e )\ R=oo d9. (45) 

Subject to the axially symmetric boundary values (fl3l) . the integrals (I4"5l) and (144]) for D = 3 
and D = 4,5,6 are evaluated as 

J 3 ax = Arfunll - (-l) m ] (46) 

I? = 12 77 4 7r 2 m (47) 

J 5 ax = 32r/ 5 7r 2 [l-(-l) m ] (48) 

J 6 ax = 5\ v 6 7c 3 m. (49) 

We now see from (T461) and (T48|) . that in odd D dimensions axially symmetric fields 
are capable of supporting zero topological charge solutions describing an even number m of 
soliton-antisoliton energy/charge concentrations, as well as unit topological charge solutions 
describing chains [H] for an odd number m. As we shall see from the numerical work in 
Section 5.2, these concentrations are located slightly off the xn axis, forming rings analogous 
to the nodes on the symmetry axis |f| found in the 3 dimensional Yang-Mills-Higgs case. We 
see by contrast from fj47l) and (|4"91) . that in even D dimensions axially symmetric fields are 
not capable of supporting zero topological charge solutions. They describe only multisoliton 
solutions of topological charges m, the concentrations of charge/energy being located on 
the xd axis. (These are the analogues of Witten's axially symmetric instantons [17].) Our 
numerical solutions in the next Section will bear out these conclusions. 

Having described candidates for zero topological charge solutions in odd dimensions, we 
proceed to explore prescriptions whereby such solutions in even D dimensions can also be 



3 Subsequent to the construction of zro charge monopole-antimonopole pairs [3 [8], such charge chains of 
monopoles and antimonopoles of unit topological charge were contructed in [9] . 
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constructed. This is possible only if less stringent symmetry than axial symmetry is imposed 
on the system, and below we describe two such distinct prescriptions in D = 4, employing in 
turn azimuthal and bi-azimuthal symmetries, and one such prescription in D = 6 employing 
tri-azimuthal symmetry. 



4.2 Azimuthal symmetry 

In the case of azimuthal symmetry in D = 4, the asymptotic boundary values to be imposed 
on the functions h(p,z,t), f(p,z,t) and g(p,z,t) defined in (133]) are 

lim h(p,z,t) = sinmi^i sinm 2 2 

lim f(p,z,t) = sinm^i cosm 2 ^ (50) 

R— >oo 

lim g(p, z, t) = cosmx^! , m 1 , m 2 = 1, 2, 3, ... 

R— >oo 

We note here that the asymptotic axially symmetric boundary values are described by 
one integer m for D > 4 and two integers (m,n) for D = 3, while those for the azimuthal 
boundary values for D = 4 are given in terms of the triple of integers (m 1; m 2 , n). 

Substituting azimuthal Ansatz (!3~T!) in (|4*2I) for D = 4, and using the analyticity re- 
quirement that h(r, t) vanishes on the t-axis, this reduces to the two dimensional angular 
integral 

If = A\n J h [h (fe 2 g dl - f 6l ge 2 ) + / (9e 2 h 6l - h dl h d2 ) + g {h e J dl - g e Je 2 )} d9 1 d9 2 , (51) 

which can readily be evaluated subject to the boundary conditions (!50|) to yield 

Jf im = Urf^mm [1 - (-l) ma ] • (52) 

This accommodates both multi-soliton (for odd m 2 ) and zero topological charge (for even 
m 2 ) solutions, labeled by the triple of integers (mi,m 2 ,n). Unfortunately the numerical 
solution the corresponding field equations involves three dimensional integration, which task 
is beyond the scope of the present work. 

4.3 Intermediate symmetries 

In both the D = 5 and D = 6 cases discussed in section 3.4 above, the asymptotic behaviours 
consistent with finite energy are both stated formally as 

lim h(r,s,t) = sinmi#i sin m 2 # 2 

_R— >oo 

lim g(r, s,t) = sinm^i cosm 2 ^ (53) 

R^oo 

lim f(r,s,t) = cosmi^i , mi, m 2 = 1, 2, 3, ... 

R^oo 

augmented by the analyticity condition h(r = 0) = 0, which is crucial in the evaluation of 
the surface intergals. In both the D = 5 and D = 6 here, these follow from three dimensional 
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volume integrals which are formally identical. Up to numerical factors, these are expressed 
as 

jmter _ ^ (n) 2 ' 3 7^ 77l 2 J 6 pup B ABC d p E A d u E B d p E° dr ds dt 

= r^V) 2 ' 3 mi ™ 2 J e pvp e ABC E A d v E B d p E c dS, , (54) 

where we have used the notation x p = (r,s,t), and, the triplet function E A , A = 1,2,3 in 
the two cases is defined in terms of the functions (h, g, f) as 

E A = ((h) 3 , g,f) and E A =((h)\g,f) , 

respectively. The nonvanishing contributions to the surface integral (s) ( |54l) come from the 
upper hemisphere. 

The values of the respective surface integrals in D = 5, 6 are calculated to be 

jinter = ^_ v 5 n 2 m2 [l_(_1)™i] ( 55 ) 

7™ ter = 1577 6 7r 3 mi [l-(-l) ma ] , (56) 

Note the roles of mi and m 2 interchanging in (I55p and (I56p . following from cancellations 
ocurring when evaluating (I54I) . 

One sees again that by relaxing axial symmetry and imposing a weaker symmetry, it 
is possible to support both multisolitons of arbitrary topological charges, and, soliton- 
antisolitons chains, with zero and nonzero topological charges in all dimensions. Unfor- 
tunately the simplest such examples result in three dimensional boundary value problems, 
which is at present technically too hard a task to perform. The situation is the same in the 
azimuthal case in D = 4 above. 



4.4 Bi- azimuthal symmetry 

Bi- azimuthal symmetry will be applied in D = 4 and D = 5, each resulting in a two and a 
three dimensional subsystems, respectively. 

L> = 4 = 2 + 2 



In this case the fields are described by the bi-azimuthal Ansatz ( l36l) . The asymptotic 
behaviours of the functions h and g in fl36|) are taken to be 

lim h = sin mijj , lim g = cos mip . (57) 

R^oo R^oo 

The topological charge in this case is 

jbi-azim = ^ 3 ; ( 27r )2 ^ ^ J £ ^ £ AB Q^A Q^B dp rf(J 

= r? 4 3!(27r)Vn 2 y (e AB E A d p E B ) ds p (58) 
15 



where we have used the notation x p = (p, a), and E A = ((h) 2 , (g) 2 ). 

Using the analyticity conditions h(ip = 0) = and g(ip = |) = leads to the vanishing 
of the line integrals on the p and the a axes, the non vanishing contribution coming from 
the infinite quarter circle contour, readily evaluated to yield 



jbi-azim = ^4 2 ^2 ^ _ ^-Qm] ? ^ 

which supports both multisolitons and zero charge soliton-antisolitons. 
£> = 5 = 2 + 2 + l 

In this case the fields are described by the bi-azimuthal Ansatz (1391) . The asymptotic 
behaviours of the functions h and g and / in (1391) are taken to be 

lim h = sin m\Q sin m-ity , lim g = smm\6 cosrri2ip , lim / = cos ni\Q . (60) 

R— >oo R^oo R— >oo 

The topological charge now reduces to a three dimensional integral in the residual coordinates 
x M = (r, s, t) 

jbi-azim = ^ g ^2 ^ ^ y £ ^ £ ABC q^A g^B q^C dr ds ^ 



= rf 5 (2nY m n 2 J e^ p e A ^ E A d v E B d p E^ dS p , (61) 

in which the triplet function E A , A = 1, 2, 3 is defined as 

E A = ((h) 2 ,(g) 2 ,f). 
The surface integral (|6"Tj) is then performed to yield 

jbi-azim = ^ 4 j ni ^ [ X _ (_!)"»!] j ( 6 2) 

describing both multisolitons and soliton antisolitons. Note that only mi, and not m 2 , 
features in (1621) . due to a cancellation in evaluating (I61"j) . 

4.5 Tri-azimuthal symmetry 

This pertains to D = 6 only. The asymptotic behaviours of the functions h, g and / in the 
Ansatz (1411) are taken to be 

lim h = smmiipismm2ip2 

lim g = sin m^i cos m^i (63) 
lim / = cosmiV'i . 

R^oo 

The topological charge integral in this case is 

jtri-azim = ^ 90 ( 2 7t) 3 n X U 2 U 3 J E pl/p E AB ° d p E A d u E B d p E° dp d(T dr 

= r/ 5 90 (2vr) 3 m n 2 n 3 ^ f> p ^ sc ~ A 0„S B «9 P ~ c , (64) 
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where we have used the notation x M = (p, a, r), and E A = ((h) 2 , (g) 2 , (f) 2 )- 

To evaluate the surface integral fl64l) we need analytic information which comes from 
finite energy conditions. While we are not displaying here the energy density functional in 
terms of the functions (h,g,f), it is nontheless easy to deduce that h(ipi = 0, ip 2 = 0) = 0. 
g(ipi = 0, ip2 = f and f(ipi = f ) = 0) = 0. These, together with continuity conditions, imply 
that the flux (IMj) out of the three quarter planes (p, a) , (a, r) and (r, p) vanishes, and hence 
the only contribution comes from the surface bounding the octant of the two-sphere with 
radius R = \J p 2 + a 2 + r 2 . 

Applying the boundary functions ( l63l) on the asymptotic octant the flux ( |64l) yields 

^ ri - azim = nm 2 n 3 Q [1 - (-1)"*]) (i [1 - (-in)' , (65) 

analogous to fl59l) . like which the topological charge vanishes when either mi or is even, 
and otherwise it is given by the product of the vortex numbers pertaining to each of the 
azimuthal symmetries imposed. 

5 Numerical constructions 

In this Section we give numerical evidence for the existence spherically symmetric, and 
axially symmetric solutions in D = 4, 5. In addition we have constructed solutions with 
bi-azimuthal symmetry in D = 4. The solutions of the corresponding D = 3 model were 
presented in [TT] to which we refer for the latter. 

Of course, the most interesting solutions from the viewpoint of understanding zero topo- 
logical charge, are the axially symmetric ones, but the spherically symmetric ones are also 
presented mainly because the equations of motion in that case allow a thorough asymptotic 
analysis underpinning the numerical work. Also the spherically symmetric solutions present 
useful starting profiles for the D = 4 bi-azimuthally symmetric multi-solitons. 

Technically, we have restricted ourselves to two dimensional numerical integration, solu- 
tions with azimuthal symmetry in D > 4 representing a difficult numerical challenge which 
we leave for future work. Also, one should notice that only one coupling constant A« is 
relevant here. For example, one may factorize Ai and, by using a suitable rescaling, one may 
set A 2 = 1 or A 3 = 1 without any loss of generality. 

To simplify the picture, in this section we shall note 9d-2 = and xd = z. Also, for 
all configurations, the total mass/energy M (which equals the total action) is computed by 
integrating the corresponding reduced energy functionals. 
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Figure 1: The scalar function Q an d weighted energy density E of two typical D = 4 and D = 5 
spherically symmetric solutions with Ai = 1 are shown as a function of the compactified radial 
coordinate R/(l + R). 



5.1 Spherically symmetric solutions 

Considering the ansatz (}26l) . the reduced one dimensional weighted energy density reads 



E = R D ^S = R - 1 [\,{Q 2 - l) 4 [Q' 2 + (D - + 2(D - 2)A 2 (Q 2 - 1) 2 |J (66) 

x (2Q' 2 + (-D — 2)g) + 6As(Z> - l)p - 2)g ( 3 Q' 2 + (D - 3)g 
which leads to the following differential equation 

2i2 fl - 1 Q' (Ax(Q 2 - l) 4 + AX 2 (D - 2)(Q 2 - l)% + 18A 3 (£> -1)(D- 2)% 

= i? D - 2 (2A 1 (Q 2 - 1) 3 (4Q' 2 + 5(D - l)g - ^ + 4A 2 ^Q(Q 2 - 1) 
+36A 3 (D - 1 ^ 2) Q 3 (2^ 2 + (-D — 3)g). 

The solutions of this equation have been constructed numerically, for a range of the parame- 
ters Aj. We follow the usual approach and, by using a standard ordinary differential equation 
solver, we evaluate the initial condition 

26 3 A 

Q(fl) - tfi - 3(A 1 + 12 y( A 2 + 9W 3 )) fl3 + 0(fl5) - forC = 4 ' < 67) 

= bR - lix^X^Xf 3 + ° {R% fM D = 5 (68) 

at R = for global tolerance 10~ 14 , adjusting for the shooting parameter b and integrating 
towards R — > oo. The behaviour of finite energy solutions for large values of R is 
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Figure 2: A three-dimensional plot of the weighted energy density E(R, 6) of a D = 5, m = 2 
axially symmetric solution with Ai = A2 = 1, A3 = 75. 

Q(R) = 1 + ce-lv/W* _ _ i_ + 0(1/jR 6 )) for D = 4? (6Q) 

2 2 

Q(J2) = 1 + ce-fv/W* - ^£-1 + 8 -±^±- 4 + 0(l/i? 6 ), for D = 5, (70) 

where c is a free parameter (the corresponding expressions for the .D = 3 model are given 
in [H]). For all cases considered, solutions with the correct asymptotics are found when the 
first derivative of the scalar function Q(R) evaluated at the origin, Q'(0) = b, takes on a 
certain value, which is a function of Aj. 

The profiles of typical D = 4, 5 solutions are presented in Figure 1 for Ai = A2 = A3 = 1. 
The weighted energy density, as given by (166]) is also exhibited (one should notice the different 
length scales of the D = 4 and D = 5 solitons). Similar to the D = 3 case, no multinode 
radial solutions were found, although we have no analytical argument for their absence. 

5.2 Axially symmetric solutions 

Scalar solitons with axial symmetry are found by taking m > 2 in the boundary conditions at 
infinity (I43p . The two-dimensional weighted energy density E(R, 0) = R D ~ l sin D_2 6 £(R, 0) 
and the set of two coupled non-linear elliptic partial differential equations satisfied by the 
functions H(R,9), G(R,9) can easily be derived by using the reduced building blocks fl30l 
and we shall not present them here. These equations are solved numerically, subject to the 
boundary conditions 

H\ R=0 = , d R G\ R=0 = 0. (71) 

at the origin and (HH at infinitjfl (we have restricted our analysis to m = 2 solutions; the 
m = 1 case corresponds to spherically symmetric configurations). Considering solutions with 

4 In the numerical algorithm we have employed a compactified radial coordinate x = R/(l +R), such that 
spatial infinity corresponds to x = 1. 
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Figure 3: The modulus of the scalar field \ip\ = V H 2 + G 2 is shown as a function of the coordinates 
r and z for a typical D = 4, m = 2 axially symmetric solution. Here r = RsinO, z = RcosO. 

parity reflection symmetry, the equations are integrated in the < 9 < n/2 region. The 
boundary conditions satisfied at the limits of the ^-interval are 

H\ e= o = dgG\o =0 = 0, dgH\ 9=7T / 2 = G\g =n / 2 = 0. (72) 

The absence of suitable starting profiles makes this problem extremely difficullEI. The nu- 
merical calculations were performed with the software package CADSOL/FIDISOL, based 
on the Newton- Raphson method [T8] . 

The numerical error for the functions is estimated to be of the order of 10~ 2 or lower for 
most of the axially symmetric configurations. 

Solutions with m = 2 of the corresponding D = 3 model were discussed in plj. In that 
case it was possible to distinguish two individual components (e.g. the modulus of the scalar 
field |0| = + 0| possesses always two distinct zeros on the z axis). 

Our D = 4, 5 results indicate that this is a generic feature of all axially symmetric 
solutions. In Figure 2 we present a three dimensional plot of the weighted energy density 
(the reduced Lagrangian) E(R, 9) of a typical D = 5, m = 2 axially symmetric solution 
as a function of the r,z (here A = A 2 = 1, A 3 = 75). The modulus of the scalar field 
\ip\ = V H 2 + G 2 of a D = 4 solution with A = A2 = 1, A3 = 8 is presented in Figure 3. 
We have found that \<p\ possesses always two zeros at ±d/2 on the z— symmetry axis, the 
positions of the nodes depending on the value of the coupling constants Aj. The total action 
of these solutions, as given by the integral of E(R, 9) increases with increasing Aj. 

Interestingly enough, the weighted energy density E(R, 9) possesses a saddle point at 
the origin, the maxima being localized at z — id/2, at a nonzero value of r, r = r . This 
feature, already present in the D = 3 case (see Figure 4 in Ref. [TT]) is enhanced for the 

5 We managed to overcome this difficulty by improving, in succesive steps, an initial guess solution con- 
structed with suitable trial functions which interpolates between the assymptotics (|7Tj) . f|43[) . 
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Figure 4: The profiles of the scalar functions h and g are shown for a typical D = 4 by-azimuthally 
symmetric solution with m = n 2 = 2, Ai = A2 = A3 = 1. 

higher dimensional configurations, in contrast with the D = 3 Yang-Mills-Higgs [8l [9] where 
the concentartions of energy are located exactly on the symmetry axis. 

Although the profiles of the axially symmetric solutions look qualitatively the same for 
D — 3, 4, 5, their physical significance is very different. For D = 4 they describe two 
distinct solitons sitting at (z = ±d/2, r = Tq), while in three and five dimensions the 
solutions represent a pair of soliton-antisoliton with zero topological charge. 

It would be interesting to construct higher m solutions, describing for an odd dimension 
soliton-anti-soliton chains, in analogy with the situation in YMH theory [9]. 

5.3 Solutions with bi-azimuthal symmetry 

To obtain D — 4, m — 1 configurations with bi-azimuthal symmetry, we employ the n — 1 
spherically symmetric solutions discussed in 5.1 above for starting profiles and increase the 
values of ni, n 2 slowly. The iterations converge, and repeating the procedure one obtains 
in this way solutions for arbitrary n. The physical values of n 1; n 2 are integers. We have 
studied solutions with 1 < rii,n 2 < 9. The weighted energy density E(R,ip) can be written 
in terms of the reduced building blocks (1371) . The two scalar functions hi^R.ip) and g(R,ip) 
satisfy the boundary conditions 

h\ R=0 = g\ R=0 = (73) 

at the origin, (!57j) at infinity, and 

h\ip =0 = d^gl^o = 0, d^h\^ =1T /2 = g\^=Tr/2 = 0. (74) 

on the p and o axes. The field equations have been solved by using the same methods 
employed in the axially symmetric case but now with much better accuracy, the typical 
numerical error being of the order of 10 -4 or smaller. 
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Figure 5: A three-dimensional plot of the weighted energy density E(p, a) of the D = 4 bi- 
azimuthally symmetric solution presented in Figure 4. 

As expected, the bi-azimuthally symmetric solutions exhibit a very different picture. A 
general feature of all m = 1 solutions with n\ = n 2 is that the weighted energy density 
E(R,ip) possesses one maximum on the ip = 7r/4 axis (corresponding to the p = o surface), 
it being possible to distinguish only one individual concentration of the action. In this 
respect the bi-azimuthally symmetric multisolitons of this model are qualitatively similar 
to those of the models featuring gauge fields, as discussed recently in [TUJ [19]. In contrast 
to the latter however, where only solutions with rii = n 2 were found, here we noticed the 
existence finite mass solutions with n\ ^ n 2 . The maximum of the weighted energy density 
moves inward with increasing rii,n 2 . 

In Figure 4 the profiles scalar functions h and g of the a typical bi-azimuthally symmetric 
solution are shown for several angles as a function of the radial coordinate R (with Ai = A2 = 
A3 = 1 in this case). A three dimensional plot of the weighted energy density of a typical 
m — 1, rti = n 2 = 2 configuration is presented in Figure 5. We have also studied the mass 
dependence of the bi-azimuthally symmetric solutions on the coupling constants Aj. From 
the numerical results, we observed some features of the solutions, without attempting to give 
an analytic explanation. The most peculiar of these is the fact that finite mass solutions 
persist for a small but finite range of negative A2 < 0. No such solutions can be justified by 
the topological lower bounds. 

When Ai or A2 are varied, the maximum of the weighted energy density moves inwards 
with increasing value of the respective coupling constant, while the opposite behaviour to 
this is found when A3 is varied. Yet another property of multi-solitons is observed. It turns 
out that the mass of a solution M(n 1; n 2 ), whose topological lower bound given by fl59l) to 
be 47r 2 n 1 n 2 , is quantitatively quite close to riin 2 M(l, 1). This means that the deviation 
of the value of M(ni,n 2 ) from its lower bound value is proportionate to the deviation of 
the value of M(l, 1) from its respective lower bound value, implying that composite solitons 
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Figure 6: The weighted energy density E is shown for a D = 4 bi-azimuthally symmetric solution 
presenting two localized elementary constituents. 

have rather low binding energies. On the other hand, in all cases studied, it turns out that 
M(ni,n 2 ) > niH2 M(l,l), albeit by a small amount quantitatively. This in turn suggests 
that the composite states are unstable against decay into the lowest mass constituents. 

Yet another surprising but not counterintuitive property observed is, that while multi- 
solitons of masses M(ni,n2) exhibit one single peak of the weighted energy density for small 
values of \n\ — n 2 \, when \ni — n 2 \ becomes large the weighted energy density develops two 
separated peaks (see Figure 6). 

More complicated bi-azimuthal solutions with m > 2 are likely to exist. These con- 
figurations would describe composite bound states, rather analogous to the monopole- 
antimonopole chains of the Yang-Mills-Higgs model [9]. Our preliminary numerical results 
already indicate the existence of zero topological charge m = 2 configurations with bi- 
azimuthal symmetry in D = 4, with rii = n 2 = 2. We found that there were no m = 2 
solutions with n\ — n 2 — 1, just as for the pure YM model whose instanton-antiinstanton 
solutions constructed in [10]. Again as in [10], the weighted energy density exhibits two 
distinct maxima on the ip = ir/4 axis. In the absence of suitable starting profiles however, 
the numerical accuracy of these solutions turned out to be much lower in this casfi 

We hope to return to a systematic discussion of these solutions, together with a general- 
isation to higher dimensions. 

6 We started the numerical process by employing a guess solution constructed with a suitable trial functions 
that contains several free parameters whose values are then tuned. However, improving the accuracy of 
these profiles has proven more difficult than in the case of axially symmetric configurations, where a similar 
approach has been employed. 
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6 Summary and Conclusions 



The overriding aim of this work is to examine the conditions that enable the construction of 
zero topological charge solutions in classical field theories, which otherwise support topolog- 
ically stable (multi-)soliton solutions with nonvanishing topological charge. We have shown 
that the conditions in question are those of (a) subjecting the system to the requisite sym- 
metry, which in practice is what is done anyway when constructing (multi-)solitons, and (b) 
by requiring special types of boundary conditions that differ essentially from those employed 
for (multi-)solitons. 

Our symmetry analysis covers dimensions D = 3, 4, 5, 6, namely both odd and even 
examples in 6 > D > 3, while the solutions constructed numerically to underpin our findings 
are limited to D < 5. The reason for this restriction is that boundary value problems in more 
than two dimensions is beyond the scope of the present work. The case D = 2 is irrelevant, 
being too small to accommodate requirements (a) and (b). The case D = 3, while it is the 
first nontrivial example, is rather special since in that case axial and azimuthal symmetries 
coincide. We have chosen to carry out this investigation in the framework of the simplest 
possible field theoretic model, irrespective of its applicability to physical problems. This 
is a symmetry breaking Goldstone type model in D Euclidean dimensions, whose energy 
density functional depends onaD component scalar field <f) a , a = 1,2, . . . , D. This choice is 
motivated by the fact that the topological charges in such models are the simplest available 
examples, being up to constant multiples the winding numbers of 4> a in R D . It should 
be emphasized however that our conclusions hold in the other classical field theories that 
support solitons in these dimensions, namely the sigma models and the non Abelian gauge 
field systems (including Higgs fields). 

Our conclusions can be summarised as follows: 

• The field a asymptotically tends to a unit vector <fi a , which depends exclusively on an- 
gular variables, the radial variable being infinite. Precisely what these angular variables 
are depends on the symmetry imposed. To analyse qualitatively distinct possibilities, 
we have found it sufficient to impose symmetries that result in residual subsystems 
of no more than three dimensions. It is superfluous to consider weaker symmetries 
resulting in four or higher number of effective degrees of freedom, since these do not 
result in qualitatively new features as far as the existence of multisolitons, and soliton- 
antisoliton chains is concerned. 

• There are two types of symmetries employed. First, spherical (rotational) symmetry 
in a N dimensional subspace of R D . For N = D — 1, this is axial symmetry resulting 
in two effective degrees of freedom. At the other extreme, N = 2, this is azimuthal 
symmetry resulting in an effective D — 2 dimensional subsystem. Accordingly, we have 
restricted to D = 3, 4 in the case of azimuthal symmetry. Intermediate values of N 
subject to this restriction are, N = 3 in D = 5 and N = 4 in D = 6. We have 
described these as intermediate symmetries. Second, we impose multi- azimuthal sym- 
metries, composed of azimuthal symmetries in pairs of coordinates. Again, subject to 
limiting our considerations to three effective degrees of freedom, these are bi-azimuthal 
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symmetry in D = 4, 5, and tri-azimuthal symmetry in D = 6. The residual system 
after symmetry imposition depends on the radial coordinate and the remaining angular 
coordinates, and the number of unknown functions is the same as the dimensionality 
of the residual subsystem. In all cases, the azimuthal angles are all integrated out and 
the remaining angular dependence is on polar angles {^} (0 < 9 < tt) and semi-polar 
angles {ipi} (0 < if)j < |-), i and / labeling the residual polar and semi-polar angles, 
respectively. 

• The asymptotic field 0" is parametrised by the residual angular variables {9i} and {ipi} 
only. The other angular variables, that include all azimuthal angles, are integrated out. 
Consistently with the requirements of finite energy and analyticity, the most general 
as (f) a is encoded by {6>j} and {ipi} is via 

{rriiOi} , {mjipi} , rrii^mj integers. (75) 

It is important to stress that the integers (m 8 , mi) in f!75|) appear only the in the 
asymptotic field a , and that they do not parametrise the field <p a everywhere. 
Throughout this text, we have reserved the letter m to these m-numbers. In con- 
trast we label the vorticity associated with each azimuthal symmetry, on which the 
field <p a everywhere depends, by the letters {n}. Thus the n-numbers which count the 
winding in each azimuthal plane are on a completely distinct footing as opposed to the 
m-numbers which serve only to select the boundary values imposed. All solutions with 
m = 1 describe topologically stable multisolitons whose topological charges are en- 
coded with the n-numbers, {n}. The topological charges of soliton-antisoliton chains 
with even m are zero, while those of odd m are nonzero, and depend on {n}. In the 
special case of axial symmetry in D > 4, when no n-number occurs, the topological 
charges in even D are labeled by a m- number. 

The numerical constructions in Section 5 underpin the above conclusions. Both 
(multi)solitons and soliton-antisoliton solutions have been constructed like for the D = 3 
case in [TT], whose results are extended to higher dimensions the present work. Various fea- 
tures found there are shared by higher dimensional axially symmetric solutions. In particular 
the profiles of the scalar functions have rather similar shapes. The numerical constructions 
in the present work is limited to 2 dimensional boundary value problem involving two func- 
tions, so that we present only axially symmetric solutions in D = 4, 5 (and D = 3 in 
and bi-azimuthally symmetric solutions in D = 4. 

The axially symmetric solutions constructed in both [TT] and in section 5.2 here are 
limited in their scope to solutions with asymptotic behaviour characterised by m— number 
equal to 1 and to 2. In the D = 3 case [11] these are multisolitons with m = 1 and higher 
n— numbers, and to soliton-antisoliton pairs with m = 2 and n— number equal to 1. In 
D = 3, like for the YMH monopoles (see e.g. [20] and [21]) the energy density of the 

7 These solutions are not absolutely stable since no first order Bogomol'nyi type are saturated in these 
models. Rather, the stability in question is a consequence of the energy respecting topological lower bound, 
like in the case of Skyrmions [l20 
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multisolitons is concentrated at the origin, and when m = 2 there occur two concentrations 
distributed symmetrically on the 3— axis. While we expect that solitons with m > 2 and 
with n = 1 would describe chains of solitons and antisolitons on the 3— axis, and when n > 3 
rings would form, like for the monopoles in the YMH model observed in [H] , this has not been 
carried out in [TT] . Here, axially symmetric solutions in D = 4, 5 are constructed in section 
5.2. The m = 2 solutions in D = 5 are similar to the m = 2 , n = 1 soliton-antisoliton 
solutions in D = 3, i.e. they describe two distinct peaks of the weighted energy density on 
the 5— axis of equal and opposite charges. But there is no n— number in D = 5 so here the 
analogy with D = 3 stops. The m = 2 axially symmetric solutions in D = 4 also describe 
two distinct particles located on the 4— axis, but unlike those of D = 3, 5 they both peaks of 
the energy density have the same topological charge. These are multisolitons, qualitatively 
similar to the axially symmetric Witten multiinstantons [17] . This illustrates that for axially 
symmetric fields with m— number higher than 1, the peaks of the weighted energy density 
are situated on the (symmetry) D— axis, such that in odd dimensions their charges have 
alternating signs, while in even dimensions all the charges have the same sign. 

There is one final property of axially symmetric solutions worth remarking on. In gauge 
field systems, namely the D = 3 YMH configurations as our only example, the zero charge 
m = 2 solutions have a positive binding energy with respect to decay into two charge— 1 
monopoles [8]. By contrast, the multisolitons and the soliton-antisoliton solutions of the 
Goldstone models have negative binding energies. Because the data available to us is not 
sharp enough, we have not displayed this quantitatively either with a plot or a table in 
section 5.2. Nevertheless the observed qualitative trend is unmistakable. 

The bi-azimuthally symmetric solutions in D = 4 constructed in section 5.3 have their 
analogue in the bi-azimuthal YM instantons given in [10]. Like in that case there is a 
n— number associated with each (of the two) azimuthal symmetries, n\ and ri2, and the 
topological charge is proportional to n\ n 2 (see (1621) ). Unlike in [TU] however, where finite 
action solutions occur only for one interger n\ = n 2 = n, here there are solutions for distinct 
rii ^ ri2. For m— number equal to 1 with n\ = ri2 = n the action density has only one 
peak which like in the YM example [10] is not situated at the origin. Rather, it peaks at 
a numerically determined distance from the origin on the ip — ^ axis, ip being the unique 
semi-polar angle. The situation is different in the (m = 1) n\ ^ n 2 case. There the action 
density breaks up into two distinct peaks on the ip — | axis, and the centres of these peaks 
move away from each other as \ri\ — ri2\ increases. Another point of contrast with the YM 
case, where the multiinstantons do form bound states, the corresponding multisolitons here 
do not form bound states. For all configurations we have studied, the energy of a the m, n 2 
multi-solitons of our model are greater than that of n\n2 1-solitons. Moreover it turns out 
that this deficit of binding energy increases with increasing m, ri2. We have also verified that 
m = 2 solutions carrying zero topological charge (see (162]) ) exist, in the n\ = n 2 = 2 case, but 
have not supplied quantitative data here. These present two distinct peaks of the weighted 
energy density like in the YM case [10] , Our numerical results here were not sufficiently 
accurate to enable us to estimate wheter the binding energy preventing the decay of this 
solution into two charge— 2 (m = 1) multisolitons is positive or negative. Likewise for the 
same reason, we did not increase n\ and n 2 to values higher than 3, to see what the analogues 
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of the rings forming in the D = 3 YM example [8] are. (Such ring like configurations were 
discovered recently in the bi-azimuthal gauge field configurations with n = 3 in |19j . implying 
their ocurrence here.) 

This completes the summary of our results. We now make some final, general comments. 
We have seen that most of the geometrical and topological properties of the multisoliton 
and soliton-antisoliton solutions in the Goldstone models studied here, are broadly similar 
both to the YMH example [8j 9| in D = 3 and to the YM example [10] in D = 4. There are 
however some notable differences, firstly that the binding energies of our multisolitons are 
negative as opposed to those of their gauge field counterparts [HI [9], [10], which are positive. 
Then there is the difference between the D = 4 bi-azimuthal Goldstone solitons, where the 
two vorticities (n— numbres) can be different, and the D = 4 YM instantons for which the two 
vorticities must be equal. More recently the SU(2) YM-dilaton system in 4 + 1 dimensions 
was analysed and the static bi-azimuthally symmetric solutions were studied in [19]. There 
too, the numerical results indicated that the two n— numbers had to be equal ni = n-i- It 
appears therefore that this restriction (n\ = n^) applies to bi-azimuthally symmetric gauge 
fields, but not to Goldstone fields. It is likely this feature may persist in multi-azimuthal 
systems too, but since this conclusion is reached only on the basis of numerics, it is beyond 
the scope of the present work. 

Based on what we have learnt about the general similarities in the different models sup- 
porting topologically nontrivial lumps studied here and in [HI [9j QUI EE] , we would speculate 
that similar analogous properties can be expected for the lump solutions in various sigma 
models, e,g. 0(D + 1) models on R D , or the corresponding Grassmannian sigma models on 
R 2iV , or indeed their gauged counterparts. One respect in which it would have been more 
appropriate to use 0(D + 1) models on R D instead of Goldstone models, featuring negative 
binding energies, is that the 0(D + 1) models would be expected to feature positive bind- 
ing energies, based on our knowledge of the 0(3 + 1) model on R 3 , namely the celebrated 
Skyrme model. Certainly the simple analysis of the topological charges and boundary con- 
ditions given in section 4 can be extended systematically and without obstacles to the sigma 
model counterparts of the scalar Goldstone fields. This was eschewed because the numerical 
constructions for the sigma models, in particular the practical task of imposing the boundary 
conditions, are very much harder. 

One last comment concerns a common feature of zero topological charge bi-azimuthal 
solutions to both gauge field systems, namely those studied in [10J [19] . and to the corre- 
spoding Goldstone model studied here. These are both solutions with m— number equal to 
2. In the former case, the numerical results indicated that the simplest such solution was 
that with n— number equal to 2, and not n = 1. Likewise in the case at hand, it turned out 
that there existed no solution for n x = n 2 = 1, the simplest solution being characterised by 
ni — 1 , ri2 — 2. It is interesting that this observation is consistent with the results of the 
numerical analysis of Krusch and Sutcliffe [22J in the context of the zero baryon charge solu- 
tions of the Skyrme model. It is very interesting also that the analytic analysis of Sadun and 
Segert |3J, which proves the existence of non-selfdual instantons (e.g. the m = 3 instantons 
in [ID]), their proof excludes topological charge 1 instantons (e.g. our m = 3 instantons with 
n = 1 whose topological charge is equal to n 2 = 1). This is a rather subtle but pervasive 
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feature, which we cannot analyse further here. 
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